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ABSTRACT 

We present an investigation into the absorber in front of the z=2.63 radio galaxy 
MRC 2025-218, using integral field spectroscopy obtained at the Very Large Telescope, 
and long slit spectroscopy obtained at the Keck 11 telescope. The properties of MRC 
2025-218 are particularly conducive to study the nature of the absorbing gas, i.e., this 
galaxy shows bright and spatially extended Lya emission, along with bright continuum 
emission from the active nucleus. 

Lya absorption is detected across ~ 40 x 30 kpc 2 , has a covering factor of ~1, and 
shows remarkably little variation in its properties across its entire spatial extent. This 
absorber is kinematically detached from the extended emission line region (EELR) . Its 
properties suggest that the absorber is outside of the EELR. We derive lower limits to 
the HI, HII and H column densities for this absorber of 3 x 10 16 , 7 x 10 17 and 2 x 10 18 
cm -2 , respectively. Moreover, the relatively bright emission from the active nucleus 
has allowed us to measure a number of metal absorption lines: CI, CII, CIV, NV, OI, 
Sill, SilV, A1II and A1III. The column density ratios are most naturally explained 
using photoionization by a hard continuum, with an ionization parameter U~0.0005- 
0.005. Shocks or photoionization by young stars cannot reproduce satisfactorily the 
measured column ratios. Using the ratio between the Sill* and Sill column densities, 
we derive a lower limit of >10 cm -3 for the electron density of the absorber. The data 
do not allow useful constraints to be placed on the metallicity of the absorber. 

We consider two possibilities for the nature of this absorber: the cosmological infall 
of gas, and an outflow driven by supernovae or the radio-jets. We find it plausible that 
the absorber around 2025-218 is in outflow. We also find good agreement between 
the observed properties of the HI absorber and the properties of the HI absorption 
expected from the cosmological infall model of Barkana & Loeb. 
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1 INTRODUCTION 

Powerful active galaxies are important probes for under- 
standing the formation and evolution of massive ellipticals 
and clusters of galaxies (e.g. Rocca-Volmerange et al. 2004; 
Kurk et al. 2000). Their prodigious luminosities allow the 
identification of large samples spanning a vast range in z 
(e.g. Rottgering et al. 1997). 



Many powerful radio galaxies at z>2 (HzRG) are em- 
bedded within giant nebulae of ionized gas which strongly 
emit both metal lines and the recombination lines of H and 
He (e.g. McCarthy et al. 1987; Villar-Martm et al. 2003). In 
addition, they often show strong absorption features in Lya 
(e.g. van Ojik et al. 1997; Wilman et al. 2004), and some- 
times in metal lines such as CIV (e.g., Binette et al. 2000; 
Jarvis et al. 2003). These absorption lines are thought to 
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be formed in cold/warm gas inbetween the emission source 
and the observer, and as such they can provide a wealth of 
information about the active nucleus, the host galaxy and 
its environment (e.g. Hamann & Ferland 1999). They are 
often spatially very extended (e.g. Rottgering et al. 1995). 
In their high spectral resolution study of the Lya profile of 
15 HzRG, van Ojik et al. (1997) noted: (i) an anticorrela- 
tion between the HI column of the absorber and the size of 
the radio source; (ii) the HI absorption features tend to be 
blueshifted relative to the line emission; and (iii) HI absorp- 
tion is more frequently detected in smaller radio sources. 
This has led various authors to suggest that the absorbers 
are part of an outflow, powered by either the radio source 
(Krause 2002; Wilman et al. 2004) or a massive starburst 
(Binette et al. 2000). In at least some HzRG, the absorbers 
are enriched in metals, although high metallicities are not 
necessarily required (Binette et al. 2006). Studying these ab- 
sorbers is clearly important for understanding the formation 
and evolution of the host galaxy, and the possible symbiosis 
between this process and the nuclear activity. 

The high-z radio galaxy MRC 2025-218 (z=2.63: Mc- 
Carthy et al. 1996) has been extensively studied in a variety 
of wavebands. Radio maps observed at 5 and 8 GHz (17 
and 30 GHz in the rest-frame; Carilli et al. 1997) show that 
MRC 2025-218 has a relatively small (-40 kpc), triple ra- 
dio source, with a jet being clearly detected on one side of 
the nucleus. Rest-frame UV and optical Hubble Space Tele- 
scope (HST hereinafter) images obtained by Pentericci et al. 
(1999; 2001) show a bright point source - the active nucleus - 
surrounded by relatively fainter, spatially extended emission 
from the host galaxy. The UV continuum has a high degree 
of linear polarization (8.3±2.3 per cent), suggesting that a 
substantial fraction of this emission is scattered nuclear light 
(Cimatti et al. 1994). In addition, UV-optical spectroscopy 
has revealed the existence of a giant nebula of warm ion- 
ized gas, with a spatial extent exceeding that of the radio 
source (Villar-Martm et al. 1999, 2007; Pentericci et al. 2001; 
Humphrey et al. 2008). This nebula emits an emission line 
spectrum characteristic of the low-density, extended emis- 
sion line regions (EELR hereinafter) of warm ionized gas 
that are commonly associated with powerful active galaxies. 
Very broad CIV and Ha emission from the broad line region 
(BLR hereinafter) is also detected (Villar-Martm et al. 1999; 
Larkin et al. 2000; Humphrey et al. 2008). Based on Spitzer 
photometry of MRC 2025-218, Seymour et al. (2007) have 
derived a upper limit to the stellar mass of M odot . 

Several properties of MRC 2025-218 are particularly 
conducive to the detection and study of intervening ab- 
sorbers. Firstly, the Lya emission from the EELR is bright 
and has a large spatial extent, allowing HI absorption to be 
traced across large spatial scales. Secondly, the brightness 
of the UV continuum emission, due to a strong contribution 
from the active nucleus itself, increases the probability of de- 
tecting lines in absorption that are usually weak or absent in 
emission in radio galaxies. In a recent paper (Villar-Martm 
et al. 2007: Paper I) we reported the detection of a strong, 
spatially extended Lya absorption feature in the spectrum 
of MRC 2025-218. Absorption from Sill, SilV, OI, CII and 
CIV have also been detected in the spectrum of this source 
(Villar-Martm et al. 1999). In this paper we present a more 
detailed analysis of the absorbing gas in front of MRC 2025- 
218, making use of integral field spectroscopy from the Very 



Large Telecope (Paper I) and a deep, long slit spectrum 
from the Keck II telescope (Humphrey et al. 2008). We also 
examine whether the absorbing gas might be related to an 
outflow of gas from the host galaxy, or the cosmological in- 
fall of gas onto the host galaxy. Throughout this paper, we 
assume a flat universe with Ho— 71 km s~ Mpc~ , Oa=0.73 
and fi m =0.27. At the redshift of this source, l"corresponds 
to 8.1 kpc. 



2 DATA 

2.1 Keck II LRIS observations 

Observations of MRC 2025-218 were made on 13 and 14 
July 1999 at the Keck II telescope using the Low Resolution 
Imaging Spectrometer (LRIS: Oke et al. 1995) in polarime- 
try mode (Goodrich et al. 1995). The LRIS detector has a 
spatial scale of 0.214"per pixel. The slit was l"wide and ori- 
ented at a position angle of 175°, i.e., along the main axis 
of the radio emission (Carilli et al. 1997). The full width 
at half of maximum (FWHM hereinafter) of the seeing disc 
was and thus the source filled the slit. A 300 line mm -1 
grating was used, providing a dispersion of 2.4 A per pixel, 
an effective spectral resolution (FWHM) of ~10 A, and a 
spectral range of ~3800-9000 A, or ~1060-2430 A in the 
rest-frame. The source was observed for a total of 5.1 hr. 
The observations were reduced using standard routines (see 
Cimatti et al. 1998; Vernet et al. 2001 for further details). 



2.2 VLT VIMOS observations 

In addition, MRC 2025-218 was observed on 2005 July 28, 
29 and 30 using the Visible MultiObject Spectrograph (VI- 
MOS: Le Fevre et al. 2003), at the UT3 unit of the Very 
Large Telescope (VLT). The integral field unit (IFU) of the 
VIMOS instrument comprises 6400 microlenses coupled to 
fibres, 1600 of which are used when in high resolution mode. 
In this mode, the IFU covers 27"x 27" on the sky and each 
fibre has a diameter of 0.67". The HRbi uc grating was used, 
resulting in an effective wavelength range of ^4150-6200 
A and an instrumental profile with a full width at half 
maximum (FWHM) of 1.7±0.2 A. The total on-object ex- 
posure time was 10 hr. During the observations, the FWHM 
of the seeing disc ranged from ~0.4"to 3.0". The FWHM of 
the final co-added data cube, measured from a bright field 
star, is 

The data were reduced using R3D (Sanchez 2006) and 
IRAF routines. The original wavelength calibration was un- 
reliable, and so the wavelength calibration was refined using 
the Keck II spectrum. In the 4650-5250 A spectral region 
we were unable to achieve an accurate wavelength calibra- 
tion, due to a lack of arc emission lines, sky lines or strong 
lines from the radio galaxy itself. See Paper I for a more 
detailed description of the observations and a discussion of 
the emission line nebula associated with MRC 2025-218. 
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Figure 1. The Keck II LRIS spectrum of MRC 2025-218, extracted from a 1.9" X 1.0"aperture. The x-axes show both the observed 
and the rest-frame wavelength scales. Identifications of the absorption lines associated with MRC 2025-218 are given. We detect two 
absorption line systems: Absorber A, with a blucshift of several hundred km s _1 relative to the narrow Hell A1640 emission line; and 
absorber B, with a blucshift of a few thousand km s _1 relative to the narrow Hell emission (see text). Unmarked absorption features 
have ambiguous identifications or are not associated with the radio galaxy. 



3 RESULTS 

3.1 Keck II LRIS spectrum 

In the LRIS spectrum, we identify two separate absorbers 
associated with MRC 2025-218. One is blueshifted by ~600 
km s _1 relative to the velocity centroid of Hell A1640, and is 
clearly detected in CI A1157, CI A1194, CI A1278, CII A1335, 
CIV AA1548,1551, NV A1239, 01 A1302, Sill A1260, Sill* 
A1265, Sill* A1309, Sill A1527, SilV A1403, A1II A1671 and 
A1III AA1855,1863 (absorber A hereinafter). Since these ab- 
sorption lines all have a similar line of sight velocity, within 
the la uncertainties, we assume that they originate from the 
same absorber. 

The other absorber (absorber B) is blueshifted by 
~2000 km s _1 relative to Hell, is detected in CI A1157, 
CI A1278, CII A1335, and SilV A1394, and has lower equiv- 
alent width (Wa) than does absorber A. Several of these 
absorption lines were previously detected by Villar-Martm 
et al. (1999) in their New Technology Telescope spectrum. 
We also detect Sill A1260 and Sill A1527 blueshifted by 3600 
and 3200 km s _1 , respectively, relative to the Hell emission; 



their blueshifts are significantly higher than those of the 
other lines from absorber B, suggesting that the Sill lines 
may come from a separate absorber. In figure 1 we show the 
LRIS spectrum with the identified absorption lines marked. 
Measurements of these absorption lines are given in tables 

1 and 2, for absorber A and B, respectively. The la errors 
associated with the equivalent width measurements were de- 
termined from the uncertainty in the height of the underly- 
ing emission. Velocity shifts are given relative to the centroid 
of the Hell emission, and equivalent widths are given in the 
rest frame of MRC 2025-218. We have not measured SilV 
A1394 from absorber A or SilV A1402.8 from absorber B, be- 
cause, due to their relative velocity shifts, these two features 
occur at very similar wavelengths. In the case of absorber 
A, we detect several excited fine structure Sill* absorption 
lines, and for this reason we assume that the CII absorp- 
tion feature contains both CII A1334.5 and the excited fine 
structure line CII* A1335.7. Also included in tables 1 and 

2 are several lines that are undetected in absorption, but 
whose upper limits are useful for constraining the column 
densities. 
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Some absorption features, although detected at the 
> 3a level of significance in the LRIS spectrum, have not 
been measured. In some cases the line identifications are 
ambiguous, i.e., when a line from absorber A and a line 
from absorber B are expected to occur at very similar wave- 
lengths. Several other absorption lines cannot be identified 
with absorber A or B, and may be formed in intervening sys- 
tems at significantly lower z. These lines are not of primary 
interest in this study, and shall notbe considered further. 

At the spectral resolution of the LRIS data, most of the 
absorption lines are unresolved in velocity space, i.e., the full 
width at half minimum is in agreement with the FWHM of 
the instrumental profile (~600 km s _1 ), within errors. For 
this reason, we consider that the FWHM of the absorber is 
likely to be <>600 km s _1 . 

The observed depth of the absorption line can be used 
to place a lower limit on the covering factor C, using the 
relation 



C> 



F 



(1) 



where Fo is the unabsorbed flux density near the wavelength 
of the absorption line, and F min is the flux density mea- 
sured at the deepest point of the absorption line profile. For 
absorber A, the deepest absorption lines (CII and CIV) re- 
quire C>0.7. In the case of absorber B, the deepest of the 
absorption lines, SilV A1394 and Sill A1527, require C>0.2. 
Since the absorption lines are not well resolved spectrally, we 
use the curve of growth (CoG hereinafter) method to derive 
column densities, or limits thereto (Spitzer 1978): when an 
absorption line is on the linear part of the CoG, the column 
density N can be derived using 

„2 



N = 



W x m e c A 
7re 2 /A 2 



(2) 



where f is the oscillator strength and Ao is the rest wave- 
length of the line; when an absorption line is saturated, 
equation 2 becomes 



N > 



W\m e c 



The implied column densities are listed in column 5 of tables 
1 and 2. While the bright emission from the active nucleus 
facilitates the detection of these absorption lines, the fact 
that this emission is spatially unresolved means that we are 
unable to determine the spatial extend of the absorption 
lines discussed above (only in the case of Lya has this been 
possible: see §3.2.1). 



3.2 VLT VIMOS spectrum 

3.2.1 Lya absorption 

The sharp edge to the blue wing in Lya, and also the fact 
that the flux drops below the continuum level, led Villar- 
Martm et al. (1999) to suggest that this line is partially ab- 
sorbed by neutral Hydrogen. Our higher spectral resolution 
VIMOS spectrum confirms this. We detect two Lya absorp- 
tion features on the blue side of the emission line (Figures 3 
and 4) with velocity shifts of ~-700 and -1700 km s^ 1 . We 
identify these two absorption features with absorber A and 
B, respectively. 

The absorption profile of Absorber A takes the form of 



a sharp absorption edge at 4402. 5±0. 3 A, with a gradual 
decrease in apparent optical depth towards shorter wave- 
lengths. At 4398.5±0.3A, there is secondary peak in the Lya 
emission, which presumably marks the wavelength at which 
absorber A becomes relatively transparent. At its deepest 
point this absorption feature is black, showing that the ab- 
sorbing gas has a covering factor of 1. 

The spatial extent of the Lya absorption from absorber 
A is strikingly large. The sharp blue edge is detected in all 
fibres which have sufficient signal to noise to detect such a 
feature (23 fibres or 15.4 arcsec 2 ). In addition, the secondary 
peak in the Lya velocity profile is detected over a similarly 
large spatial extent (16 fibres or 10.7 arcsec 2 ). The wave- 
length of the blue edge, and that of the secondary emission 
peak, are remarkably constant from fibre to fibre, varying by 
less than 100 km s _1 . The combined spatial extent of these 
two features allows us to set a lower limit of 4.7" x 3.4" 
(~40x30 kpc 2 ) for the projected size of this Lya absorber. 

Due to the complexity of the Lya velocity profile, it is 
necessary to fit the line using model velocity profiles. Our 
model profiles consist of a Gaussian shaped Lya emission 
component, a linear continuum, and an absorption compo- 
nent which has a sharp absorption edge and a Voigt-like 
decay in column density towards bluer wavelengths. The 
covering factor of the absorber, as seen by the Lya neb- 
ula and the continuum, was set to 1. Such a large covering 
factor is justified by the fact that the absorption feature is 
essentially black at the absorption edge. All other parame- 
ters were allowed to vary freely, within physically reasonable 
limits (e.g., emission and absorber wavelengths: 4360-4440A; 
emission FWHM: 100-3000 km s _1 ; absorber velocity width: 
10-3000 km s _1 ; HI column: 0-10 22 cm -2 ). We have defined 
the width of the absorber as the distance in velocity space 
from the absorption edge to where the HI column has half 
its peak value. 

Since we are also interested in the spatial properties of 
this absorption feature, we have fitted the Lya profile in four 
apertures: (i) the 'nuclear' aperture, comprising 3x3 fibres 
(i.e. 2" x 2") centred on the fibre wherein the very broad CIV 
emission and the continuum emission are brightest, which we 
assume marks the projected position of the active galactic 
nucleus (AGN); (ii) the 'Northern' aperture, defined as 3x3 
fibres centred on the fibre in which Lya is brightest; (iii) 
the 'Southern' aperture, which comprises 3x3 fibres with a 
centre that is 2 fibres South of the continuum peak; and (iv) 
the 'spatially integrated' aperture, which is the sum of all 21 
fibres contained within the above apertures. The apertures 
were sized in order to maximise the signal-to-noise ratio in 
the line profile. Note that the nuclear aperture is not inde- 
pendent of the Northern and Southern apertures. In figure 2 
we show the continuum-subtracted Lya image of 2025-218, 
reconstructed from the VIMOS data cube (see paper for 
further details of this image), with extraction apertures in- 
dicated. We have assumed that the bright continuum point- 
source (Pentericci et al. 1999; 2001) and the very broad CIV 
emission (-9000 km s _1 : Villar-Martm et al. 1999; figure 3 
of this paper) mark the position of the active galactic nu- 
cleus (AGN hereinafter). 

The fits to Lya are shown in Fig 3, and the best-fitting 
parameters are detailed in Table 3. However, we must em- 
phasise that these fits are not unique: since the absorption 
feature is saturated (i.e., black) at some velocities, the true 
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Table 1. Absorber A: measurements and limits from the LRIS 
spectrum. Columns: (1) atom/ion; (2) rest-frame wavelength (A); 
(3) velocity shift (km s _1 ) relative to the central wavelength of 
the Hell A1640 emission line; negative values indicate a blueward 
shift; lcr uncertainties are typically 100-200 km s _1 ; (4) rest- frame 
equivalent width (A); (5) column density of the atom/ion. 



Atom/ion 




Av 


Wrest 


N 




(A) 


(km s- 1 ) 


(A) 


(10 14 cm- 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


CI 


1157.2 


-670 


2.4±0.8 


>2.5 


CI 


1193.5 


-880 


1.9±1.5 


>5.5 


CI 


1277.5 


-790 


0.6±0.3 


>1.9 


CI 


1561.1 




<0.4 


<3.0 


CI 


1657.2 




<0.4 


<1.2 


CII,CII* 


1335.3 


-650 


2.0±0.2 


>8.9 


CIV 


1549.5 


-560 


4.7i0.9 


>6.3 


NV 


1238.8 


-460 


0.8±0.3 


>2.2 


OI 


1302.2 


-500 


0.8±0.3 


>5.7 


Sil 


1845.5 




<0.3 


<0.37 


Sill 


1260.4 


-620 


0.8±0.2 


>0.43 


Sill* 


1264.7 


-580 


0.3±0.2 


>0.08 


Sill* 


1309.3 


-640 


0.5±0.2 


>2.1 


Sill 


1526.7 


-590 


0.7±0.1 


>2.2 


Sill 


1808.0 




<0.15 


<24 


Sill* 


1817.2 




<0.36 


<54 


SilV 


1402.8 


-550 


1.2±0.2 


>2.3 


A1II 


1670.8 


-630 


0.8±0.1 


>0.16 


A1III 


1854.7 


-550 


0.7±0.1 


0.47±0.06 


Aim 


1862.8 


-670 


0.3±0.1 


0.35±0.12 



Table 2. Absorber B: measurements and limits from the LRIS 
spectrum. Columns: (1) atom/ion; (2) rest-frame wavelength (A); 
(3) velocity shift (km s~ 1 ) relative to the central wavelength of 
the Hell A1640 emission line; negative values indicate a blueward 
shift; lcr uncertainties are typically 100-200 kms ; (4) rest-frame 
equivalent width (A); (5) column density of the atom/ion. The 
large velocity shift between Sill and the other lines suggests that 
the Sill lines may come from a separate absorber. 



Atom/ion Ao Av W re st N 

(A) (km s- 1 ) (A) (10 14 cm- 2 ) 



(1) 


(2) 


(3) 


(4) 


(5) 


CI 


1157.2 


-2500 


1.6±0.7 


>1.4 


CI 


1277.5 


-2200 


0.6±0.2 


>2.5 


CII 


1334.5 


-2200 


0.5±0.2 


>1.5 


Sil 


1845.5 




<0.36 


<0.65 


Sill 


1260.4 


-3600 


0.7±0.3 


>0.3 


Sill 


1526.7 


-3200 


0.46±0.05 


>1.5 


SilV 


1393.8 


-2500 


0.8±0.1 


>0.80 



HI column densities could be substantially higher than these 
values. Therefore, we adopt a lower limit of 10 14 ' 8 cm~ 2 for 
the HI column of absorber A. This takes into account the 
formal uncertainty from the fitting, and represents the min- 
imum possible HI column density for the absorber. In addi- 
tion, the fact that the Lya absorption feature is saturated 
means the HI column density could vary significantly be- 
tween apertures with no obvious difference in the observed 
absorption feature. In our best fit model, the absorption edge 
is blueshifted by 380 km s _1 relative to the Hell emission. 
The FWHM of the unabsorbed Lya profiles (1200-1600 km 



s _1 ), are relatively large in comparison with that of Hell 
(~500-1100 km s _1 ). For the velocity width of the absorber 
itself, we obtain values in the range 250-350 km s _1 . 

Lya absorption from absorber B is also detected (Fig- 
ures 3 and 5), but it is not clear whether it is spatially 
extended. Given the depth of this absorption feature, we 
estimate its covering factor C to be J>0.7. It has an equiv- 
alent width of 1.9±0.3A, implying a HI column density of 
>3.5xl0 14 cm- 2 . 

Absorbers A appears to be kinematically distinct from 
the line emitting gas that comprises the EELR. While the 
line emitting gas shows a clear velocity gradient across the 
Lya nebula, the velocities at which absorber A occurs does 
not show any significant spatial variation. In addition, the 
fact that absorbers A and B both have C~l implies that 
the absorbers are very probably outside of the Lya emitting 
region. Since Lya emission is detected at distances of up to 
25 kpc from the AGN (Villar-Martm et al. 2007), we will 
assume that absorbers A and B are at distances of >25 kpc 
from the AGN. 

3.2.2 CIV absorption 

The CIV velocity profile (see Figs 4 and 5) also shows a sharp 
blue edge, which occurs at 5527 A, and which we identify 
as CIV A1550.8 from absorber A. This absorption feature 
is not black at its deepest point, and its depth implies r~l 
and a CIV column of ^3xl0 15 cm~ 2 . At ~5618 A the 
velocity profile shows a peak, and blueward of this there is a 
second absorption feature, wherein the flux reached zero. We 
identify this feature as CIV A1548.2 absorption by absorber 
A. The relative wavelengths of these two absorption features 
are in good agreement with the theoretical separation of 
the CIV doublet. The apparent depth of the CIV A1548.2 
absorption line requires tJ;2 and a CIV column of Jj3xl0 14 
cm~ 2 . Thus, independently of model profile fitting, we can 
state that the column density of CIV is in the range 10 14 ' 5 - 
lO 15 ' 5 cm' 2 . 

The velocity profile of CIV shows an additional peak 
blueshifted by ~1500 km s _1 relative to Hell. Blueward of 
this peak, at A =5553-5587 A (corresponding to ~-1500 
to -3600 km s _1 relative to Hell), the flux dips below the 
continuum level. This suggests that some of the flux within 
this wavelength range has been absorbed, possibly by CIV 
associated with absorber B. Unfortunately, the relatively low 
signal to noise ratio, along with the presence of a strong 
sky-subtraction residual at ~5577 A, precludes a detailed 
analysis of the possible CIV absorption from absorber B. 
The signal to noise of our data is insufficient to determine 
whether the CIV absorption features are spatially extended. 

To place further constraints on the CIV column of ab- 
sorber A, we have fitted model velocity profiles to the ob- 
served profile of the CIV doublet. We restrict this analysis 
to the nuclear aperture, where the signal to noise is rela- 
tively high and where the CIV absorption feature is clearly 
detected. In our fits, we have used two Gaussians for the 
CIV emission; one to represent the narrow CIV emission 
doublet from the EELR, and the other one to represent the 
CIV emission from the BLR* that is very clearly seen in our 

* Treating the emission components as single lines instead of dou- 
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3.2.3 Other absorption lines in the VIMOS spectrum 



-3 -6 -4 -2 2 4 6 8 

Offset in RA (") 

Figure 2. Continuum-subtracted Lya image of MRC 2025-218, 
reconstructed from the VIMOS data cube (Paper I). The nu- 
clear, northern and southern apertures (see text) are marked with 
boxes. North is toward the top, East is to the left. The spatial 
zero of the axes corresponds to the position of the active nucleus. 



VIMOS spectrum (Fig 3: FWHM~9000 km s _1 ; see also the 
spectrum of Villar-Martm et al. 1999). The FWHM of the 
EELR emission was set to the FWHM we have measured 
for Hell in the nuclear aperture (1100 km s _1 ). The other 
parameters of the line emission components were allowed to 
vary freely. 

We fitted the absorption component of absorber A as 
a doublet, i.e. CIV A1548 and CIV A1551, using a pair of 
profiles with a sharp absorption edge and a Voigt-like decay 
towards shorter wavelengths. The two absorption profiles 
were forced to have the same velocity shift, velocity width 
and New- In order to reduce the number of free parameters, 
the velocity width of the CIV absorber was fixed to the value 
obtained from the fit to Lya in the nuclear aperture. Nciv 
and the velocity shift were allowed to vary freely. We set 
the weighting of the data between 5553 A and 5587 A to 
zero, because this region contains a strong sky subtraction 
residual, as well as possible CIV absorption from absorber 
B and possible Sill* A1533 absorption from A. 

The best fit to the CIV velocity profile is shown in Fig 
4. Table 4 lists the parameters of the best fitting model. We 
find that N C /v = 10 15 ' 3±0 ' 3 cm" 2 . The velocity shift of the 
CIV A1551 absorption edge (-150 km s _1 ) is in agreement 
with that of the Lya edge, after taking into account the ~200 
km s _1 uncertainty in the relative wavelength calibration 
between the Lya and CIV regions of the spectrum. Due 
to the complexity of the observed CIV velocity profile, our 
best fit parameters for the emission components may not be 
unique. However, after repeating the fits with the FWHM 
and wavelength of the emission components set to a variety 
of different values, we found that our best fit parameters to 
the absorption doublet are not significantly affected by the 
precise profile of the CIV emission: this is due to the high 
optical depth of the CIV absorption lines. 



blcts has no significant impact on our model velocity profiles, be- 
cause the FWHM of both emission components are substantially 
larger than the theoretical doublet separation. 



In addition to Lya and CIV, we have detected several other 
absorption lines in our VIMOS spectra of MRC 2025-218. 
From absorber A, we detect absorption due to NV A1239, 
Sill A1260, CII+CII* A1335, SilV A1403 and Sill A1526. 
From absorber B we detect NV A1239, CI A1278 and SilV 
A1394. We also detect Sill A1260 and Sill A1526 absorp- 
tion blushifted by 3600 and 3200 km s" 1 , respectively, rel- 
ative to the Hell emissln a recent paper (Villar-Martm et 
al. 2007: Paper I) we reported the detection of a strong, 
spatially extended Lya absorption feature in the spectrum 
of MRC 2025-218. Absorption from Sill, SilV, OI, CII and 
CIV have also been detected in the spectrum of this source 
(Villar-Martm et al. 1999). In this paper we present a more 
detailed analysis of the absorbing gas in front of MRC 2025- 
218, making use of integral field spectroscopy from the Very 
Large Telecope (Paper I) and a deep, long slit spectrum 
from the Keck II telescope (Humphrey et al. 2008). We also 
examine whether the absorbing gas might be related to an 
outflow of gas from the host galaxy, or the cosmological in- 
fall of gas onto the host galaxy. Throughout this paper, we 
assume a flat universe with Ho— 71 km s _1 Mpc -1 , £7a=0.73 
and i! m =0.27. At the redshift of this source, l"corresponds 
to 8.1 kpc. 

ion. As discussed earlier, it is not clear whether these 
lines are from absorber B or from a separate one. In Fig- 
ures 5 and 6 we show two spectral regions of the VIMOS 
spectrum, with absorption line identifications marked, along 
with corresponding regions of the LRIS spectrum for com- 
parison. (We do not show the 4650-5390 A region of the 
VIMOS data because its wavelength calibration is unreli- 
able: see §2.2.) We simply measure the equivalent width of 
each absorption line, and obtain column densities, or limits 
thereto, using the CoG method (§3.1). Tables 5 and 6 list the 
absorption line properties and the implied column densities. 
As we discussed in §3.1, we are unable to determine the spa- 
tial extent of these absorption lines because the background 
source, the active nucleus, is not spatialy resolved. 

3.3 Final column densities of A and B 

For some species we have several measurements, or limits, 
of column density; where this is the case, we have reduced 
multiple measurements and limits to a 'final' column density 
as follows. (Unless otherwise stated, the procedures below 
were applied to both absorber A and absorber B.) 

• For both NV and Sill, the possible range in column 
density is constrained by lower and upper limits. In the case 
of NV, our adopted column density corresponds to the mid- 
point between the upper and the lower limit. For Sill, we 
use the midpoint between the most stringent upper limit and 
the most stringent lower limit. In both cases, we adopt a la 
uncertainty of half the difference between the two limits. 

• In the cases of CI and Sill*, the most stringent upper 
and lower limits are in fair agreement (i.e., they differ by less 
than a factor of 2), though the upper limit is actually below 
the lower limit. As before, we adopt the midpoint between 
these two limits, and take the la uncertainty to be half the 
difference between the two limits. 

• For species with either multiple upper limits or multiple 
lower limits, we use the most stringent of these limits. 
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Figure 3. Fits to the Lya profile of MRC 2025-218: (a) velocity profile in the spatially integrated aperture; (b) velocity profile in the 
nuclear aperture; (c) velocity Northern aperture; (d) velocity profile in the Southern aperture. 



• In the case of AMI, we simply take the mean of our two 
column density measurements. 



Nhi k, 3 x 10 16 . In the case of absorber B, we only have a 
lower limit to Nci, thus we cannot apply this method. 



In principle, comparing the column densities of the 
metallic species against that of HI can provide constraints 
on the metallicity of the absorbing gas. We assume that the 
column density of neutral Carbon relative to that of neu- 
tral Hydrogen is roughly equal to the Carbon to Hydrogen 
abundance ratio, i.e., Nci/Nhi ~ Nc/Nh- In the case of 
absorber A, the value of Nci and the Nhi value from our 
fits to the Lya profile imply a C/H abundance ratio that 
is Ss600 times the solar ratio of Anders & Grevesse (1989). 
Clearly, C/H abundance ratios of several hundred times the 
solar value would be unrealistic in the context of gas at a 
distance of several tens of kpc from the galaxy nucleus; as 
we suggested in §3.2.1. 

Conversely, by adopting a value for the C /H abundance 
ratio, we can then use our measurement of Nci to refine 
our lower limit on Nhi from §3.2.1. We assume conserva- 
tively that this abundance ratio is less than ten times its 
solar value. We believe this assumption to be physically rea- 
sonable in the context of spatially very extended gas out- 
side of the galaxy nucleus. As before, we also assume that 
Nci/Nhi ~ Nc/Nh - The measured CI column then implies 



We can apply a similar method to obtain a lower 
limit for the column density of ionized Hydrogen, using 
our measurement of Nsiv- It is reasonable to assume that 
in a region of warm ionized gas, the ionization fraction 
SilV/Si does not exceed the Hydrogen ionization fraction 
HII/H, i.e., Nsnv /Nsi < Nhii/Nh- If we also assume that 
Nsi/Nn <10 times its solar value, the SilV column density 
measurements then imply Nhii > 7.3xl0 17 and > 2.6xl0 17 
cm' 2 for absorber A and B, respectively. 

We can also set a lower limit to the total column of 
Hydrogen. For a given metal, we sum the column densities 
of the different ionization stages, and we then calculate the 
Hydrogen column density that is implied by the summed 
column density. As before, we assume the metal abundances 
are no higher than 10 times their solar values. In the case 
of absorber A, the columns of Silicon and Aluminium both 
imply Nh ^ 2 x 10 18 cm -2 , which we adopt as our final 
value of the total Hydrogen column of this absorber. (The 
summed column of Carbon implies Nh il 3 x 10 17 cm -2 .) 
For absorber B, the summed Silicon column implies Nh ^ 
7 x 10 18 cm -2 . In Table 7 we list our final column densities. 
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3.4 Mass and density 

In order to obtain crude estimates for the masses absorbers 
A and B, we assume that both absorbers are spherically sym- 
metric shells, and hence M H = 4:irr 2 mHN H ■ Also assuming 
that absorber A is at a distance of >25 kpc from the AGN 
(§3.2.1), and adopting N H k, 2 x 10 18 cm" 2 (§3.3), we esti- 
mate its mass to be > 1O 8 M0. In the case of absorber B, we 
also adopt a distance of >25 kpc from the AGN, and hence 
we obtain a mass of > 4 x 1O 8 M . However, a different set of 
assumptions could result in substantially different masses. 

The ratio of the Sill and Sill* column densities (i.e., 
the population ratio between the 3s 2 3p 2 Pi/ 2 and 3s 2 3p 2 P 3 / 2 
levels of Sill) is sensitive to electron density n e , and hence 
can be used to obtain an estimate of the density of the ab- 
sorbing gas. For absorber A, the SiII*/SiII column ratio is 
>0.04 (see table 7). Comparing this value against calcula- 
tions of column ratio as a function of n e (Keenan et al. 
1985), we find that n e kW cm -3 . We will assume that the 
number density of Hydrogen hh has a similar value to n e . 
In the case of absorber B, our upper limit to the SiII*/SiII 
column ratio implies n e <1000 cm -3 . Although excitation 
by 35/xm IR photons can also populate the 3s 2 3p 2 P 3 / 2 ex- 
cited level, it seems unlikely that there would be enough IR 
radiation to populate the excited level in this way. For in- 
stance, if the IR radiation from the AGN and host galaxy 
were responsible for populating the 3s 2 3p 2 P 3 / 2 Sill* level, 
then a specific luminosity at 35/im of ~ 10 39 erg s _1 Hz -1 
would be required. However, the specific IR luminosities of 
powerful active galaxies at z~2-3 are more than 4 orders of 
magnitude lower than this (e.g. Seymour et al. 2007). 

3.5 Ionization of absorber A 

Several previous studies have considered the possible source 
of ionization of the absorbers around high-z radio galax- 
ies. Binette et al. (2006) performed detailed photoionization 
calculations for the spatially extended absorbers around two 
other high-z radio galaxies, with the goal of explaining the 
measured Nciv /Nhi ratio. They found that hot stars, or a 
metagalactic background radiation (MBR hereinafter) field 
in which stars dominate over AGN, are equally successful 
for reproducing this column ratio. We now consider various 
sources of ionization for absorber A. 

Young stars. The measured Nnv /Nciv ratio of 1.8l 3 ' 2 re- 
quires a hard ionizing continuum, and hence we reject young 
stars as the main source of ionization. 

Metagalactic background. Adopting the intensity J v = 
10~ 21 erg cm~ 2 s" 1 Hz -1 sr _1 inferred by Cooke et al. 
(1997), and using the energy distribution of Fardal, Giroux 
& Schull (1998), we estimate that the MBR has an ionizing 
photon flux of ~ 2 x 10 6 cm~ 2 s _1 . This is roughly in agree- 
ment with the flux needed to explain the HII column of the 
absorber (a rec n e NHii iJ 10 6 cm~ 2 s _1 ). Thus, the ioniza- 
tion of absorber A can be explained with photoionization by 
the MBR. 

Active nucleus. The clear detection of an unresolved con- 
tinuum source in rest-frame UV and optical images of MRC 
2025-218 (e.g. Pentericci et al. 1999; 2001; also Paper I) in- 
dicates that the active nucleus is not completely obscured 
along our line of sight. Might the radiation field of the ac- 
tive nucleus be responsible for the ionization of absorber A? 



We extrapolate the ionizing photon flux from the observed 
continuum flux at 1260 A, adopting the spectral energy dis- 
tribution of Korista et al. (1997), with the exponential cut 
off set to kT —S6 eV, and assuming that the continuum is 
reddened with A„ >1.4 mag (Larkin et al. 2000). We es- 
timate the ionizing photon flux r=25 kpc from the active 
nucleus to be < 10 8 s _1 cm~ 2 , compared to the flux of 10 6 
cm~ 2 s _1 required to explain the HII column of absorber A. 
Thus, the radiation field of the active nucleus could also be 
responsible for the ionization of absorber A. 
Shocks. In explaining the ionized gas associated with active 
galaxies, shocks are often considered as an alternative to 
photoionization (e.g. Dopita & Sutherland 1995). Here, we 
compare the column density ratios measured in absorber A 
against those predicted by the fast shock models of Dopita 
& Sutherland (1996: DS96 hereinafter). We find that the 
pure shock models (i.e., ignoring the photoionized precursor 
region) are unable to explain the column ratios measured in 
absorber A: these models predict a much lower ionization 
state than is implied by our measurements, i.e., in absorber 
A we measure Nsuv/Nsui > 0.1, but in the DS96 models 
the N suv /N sui ratio is at least 100 times lower. 
Shock continuum. Although pure shocks are unable to 
explain the column density ratios of absorber A, it is in- 
teresting to consider whether the strong UV radiation field 
emitted by hot, post shock cooling plasma (e.g. Daltabuit 
& Cox 1972; Dopita, Binette & Schwartz 1982) might be 
able to explain the ionization of absorber A. To explain the 
HII column of the absorber, the shock would need to have 
a surface area of k 1.5 x 10 43 cm 2 (or k 1.6 kpc 2 ), and 
contain a gas mass of <; 1O 7 M (see DS96). To put this 
into perspective, we estimate that the EELR contains ~ 10 9 
Mq of warm ionized gas, based on the Lya lumnosity of the 
EELR (e.g. Villar-Martm et al. 2003) . Is photoionization by 
shocked gas able to explain the column ratios of absorber A? 
To explore this, we have calculated photoionization models 
using the multipurpose code MAPPINGS Ic (Binette, Dopita 
& Tuohy 1985; Ferruit et al. 1997), consisting of an iso- 
baric, plane parallel slab of gas onto which an ionizing con- 
tinuum impinges. We adopt a spectral energy distribution of 
F v oc w -1 ' 5 as an approximation of the ionizing continuum 
emitted by the shocked gas [see Fig 3 of DS96; our results 
would not be significantly affected if we use a slightly softer 
(F v oc i>~ 2 '°) or slightly harder (F v oc v~ 10 ) energy distribu- 
tion] . In order to satisfy simultaneously the observed column 
ratios Nciv/Ncii < 0.4 and Nsuv /Nsui > 0.11, an ion- 
ization parameter of U = Q/(4-Kd 2 n H c) = 0.0005 — 0.005 
would be required, implying Nhii = 7.5 x 10 19 to 7.5 x 10 19 
cm~ 2 and an ionizing flux of 5 x 10 9 to 5 x 10 10 cm~ 2 s~\ 
consistent with the lower limits derived earlier in this paper. 

In summary, there are three plausible explanations for 
the ionization of absorber A: the metagalactic background 
radiation; the radiation field of the active nucleus; or con- 
tinuum emitted by shocked gas. 



4 ORIGIN OF ABSORBER A 

Understanding the precise nature of the absorbing gas 
around high-z galaxies requires knowledge of the systemic 
velocity of the host galaxy. Since we do not detect photo- 
spheric lines in the spectrum of MRC 2025-218, we do not 
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Table 3. Results from fitting the Lyc* profile. Columns: (1) name of aperture; (2) central wavelength of emission Gaussian (A); (3) 
FWHM of emission Gaussian (km s~ 1 ); (4) wavelength of the peak column density of the absorption component (A); (5) column density 
of the absorption component in logarithm (cm -2 ); (6) Velocity width of the absorption component, defined as the distance in velocity 
space from the absorption edge to the velocity where the HI column reaches half its peak value (km s _1 ); (7) Velocity shift, relative to 
the central wavelength of the Hell A1640 emission line (km s — 1 ), of the peak optical depth of the Lyo absorber; the la errors take into 
account fitting uncertainties and the uncertainty in the wavelength calibration; the fitting uncertainties alone are 20-30 km s _1 . 



Aperture 


A e m. 


FWHMem 


Kb 


Log N HI 


Width 


Av a6 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Integrated 


4405±1 


1600±200 


4402.1±0.3 


15.1±0.3 


300±60 


-380±160 


Nucleus 


4404±1 


1600±200 


4402.2±0.3 


15.1±0.3 


350±60 


-380±160 


North 


4409±1 


1500±200 


4402.1±0.4 


15.1±0.3 


250±60 


-380±160 


South 


4402±1 


1200±200 


4402.0±0.3 


15.2±0.3 


300±50 


-390±160 



Table 4. Results from fitting the CIV profile in the nuclear aperture. Columns: (1) central wavelength of BLR emission Gaussian (A); 
(2) FWHM of BLR emission Gaussian (km s _1 ); (3) central wavelength of EELR emission Gaussian (A); (4) FWHM of EELR emission 
Gaussian (km s _1 ); (5) central wavelength of the absorption component (A); (6) column density of the absorption component (cm -2 ) 
in Logarithm; (7) Velocity width of the absorption component (km s _1 ; constrained from the Lya fit); (8) Velocity shift, relative to the 
central wavelength of the Hell A1640 emission line (km s — 1 ), of the peak optical depth of the CIV absorber; the uncertainty takes into 
account the uncertainty in the fit, and also the uncertainty in the wavelength calibration; the fitting uncertainty is 32 km s . 



A si fl FWHM BLfl Xeelr FWHMeeih Kb Log TSS C iv Width Av a6 

(1) (2) (3) (4) (5) (6) (7) (8) 



5606±5 9000±800 5614±2 1100 5615.3±0.6 15.3±0.2 350 -150±130 
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Figure 4. Fit to the CIV profile of MRC 2025-218 in the nuclear 
aperture. The weighting of the data between 5520 A and 5590 
A was reduced, because this region is strongly affected by residu- 
als from the sky-subtraction (this includes a large spike at ~5580 

A). 

know the systemic velocity of the host galaxy, and thus we 
are unable to distinguish whether absorber A is in infall or 
outflow. Therefore, both possibilities will be considered. As 
is usually the case with absorbers in front radio galaxies, we 
are unable to rule out that this absorber is not associated 
with the MRC 2025-218; however, in the interest of brevity 
we will not give further consideration to this particular pos- 
sibility. 

4.1 Outflow 

Much evidence has accumulated in support of the idea that 
part of the ionized gas of some high-z radio galaxies is in 
outflow. The wedge-shaped Lya-emitting structure oriented 



4300 4400 4500 4600 




1180 1200 1220 1240 1260 1280 

Wavelength (A) 

4300 4400 4500 4600 




1180 1200 1220 1240 1260 1280 

Wavelength (A) 

Figure 5. Absorption lines detected in the blue part of the VI- 
MOS spectrum (upper panel); the spectrum was extracted from 
the fibre associated with the nucleus, to maximise the signal to 
noise. For comparison, we also show the corresponding spectral 
region of the LRIS spectrum (lower panel) . 



perpendicularly to the radio axis in TN J1338-1942 (z=4.1; 
Zirm et al. 2005) and the extreme emission line kinemat- 
ics in the EELR of 1138-262 (z=2.2) are particularly spec- 
tacular examples of this phenomenon. The preponderance 
of blue asymmetries in the emission lines provides further 
evidence for the outflow of ionized gas from high-z radio 
galaxies (Jarvis et al. 2003; Humphrey et al. 2006). 

A number of previous investigator have proposed sce- 
narios in which the absorbers around high-z radio galaxies 
are also in outflow (see e.g. Krause 2002; Jarvis et al. 2003; 
Wilman et al. 2004). In this section, we will examine whether 
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Figure 6. Absorption lines detected in the red part of the VI- 
MOS spectrum (upper panel); the spectrum is from the Nuclear 
aperture. For comparison, we also show the corresponding spec- 
tral region of the LRIS spectrum (lower panel). 



Table 6. Absorber B: measurements and limits for other lines 
in the VIMOS spectrum. Columns: (1) atom/ion; (2) rest-frame 
wavelength (A); (3) velocity shift (km s — 1 ) relative to the cen- 
tral wavelength of the Hell A1640 emission line; negative values 
indicate a blueward shift; la uncertainties are 100-200 km S ; 
SilV falls in the region of the spectrum where the wavelength 
calibration is unreliable, and for this reason we have not calcu- 
lated its velocity shift; (4) rest-frame equivalent width (A); (5) 
column density of the atom/ion. The large velocity shift between 
Sill and the other lines suggests that the Sill lines may come from 
a separate absorber. 



Atom / ion 


An 


Av 


w rest 


N 




(A) 


(km s _1 ) 


(A) 


(10 14 cm" 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


HI 


1215.7 


-1700 


1.9±0.3 


>3.5 


CI 


1277.5 


-2200 


0.9±0.1 


>5.0 


NV 


1238.8 


-2300 


0.8±0.3 


>2.3 


Sill 


1260.4 


-3700 


1.1±0.1 


>0.71 


Sill* 


1264.7 




<0.9 


<0.70 


Sill 


1526.7 


-3200 


0.6±0.3 


>1.1 


SilV 


1393.8 




1.0±0.2 


>0.92 



Table 5. Absorber A: measurements and limits for other lines 
in the VIMOS spectrum. Columns: (1) ion; (2) rest-frame wave- 
length (A); (3) velocity shift (km s -1 ) relative to the central 
wavelength of the Hell A1640 emission line; negative values indi- 
cate a blueward shift; 1<t uncertainties are 100-200 km s - 1 ; the 
CII,CII* and SilV lines fall in the region of the spectrum where 
the wavelength calibration is inreliable - for these lines reason 
we have not calculated velocity shifts; (4) rest-frame equivalent 
width (A); (5) column density of the ion. 



Ion A() Av Wrest N 

(A) (km s" 1 ) (A) (10 14 cm~ 2 ) 

(1) (2) (3) (4) (5) 



CII.CII* 


1335.3 




1.9±0.2 


>8.4 


NV 


1238.8 


-570 


0.8±0.4 


>2.2 


NV 


1242.8 




<0.4 


<4.9 


Sill 


1260.4 


-690 


1.2±0.3 


>0.28 


Sill* 


1264.7 




<1.0 


<1.2 


Sill 


1526.7 


-460 


0.7±0.1 


>2.2 


SilV 


1402.8 




1.4±0.2 


>2.7 



it is feasible that absorber A is in outflow from the host 
galaxy. 

The kinematics of the absorbing gas appear consistent 
with those of starburst driven outflows observed in absorp- 
tion in star forming galaxies. Specifically, we observe a sharp 
absorption edge in the blue wing of Lya and CIV, where the 
flux goes to zero, with Lya also showing a gradual decrease 
in the apparent optical depth towards shorter wavelengths 
(cf. Conti et al. 1996; Pettini et al. 2000). 

In the simplest scenario, supernova explosions associ- 
ated with a massive starburst produce an expanding bubble 
of hot gas. The superbubble itself is expected to contain hot 
gas (T~ 10 6 K) which, due to its high level of ionization, 
would not be observable in the UV absorption lines consid- 
ered herein. At some point in its expansion, the gas at the 
boundary of the bubble comes into pressure equilibrium with 
the ambient intergalactic medium (IGM), and cools very ef- 
ficiently thereafter. This cooled gas would be detectable in 
various UV absorption lines. In the case of MRC 2025-218, 



Table 7. Physical properties of absorbers A and B. The first 15 
lines list the column densities, which are given in units of 10 14 
cm -2 . The final four lines give the distance of the absorber from 
the AGN, the thickness of the absorber, its mass and its density. 



Atom/ion 


A 


B 


HI 


^300 


^3.5 


HII 


^7300 


^2600 


H 


^20000 


^70000 


CI 


1.4±0.3 


>5 


CII 


>9.9 


>1.5 


CIV 


2 0+ 20 
z - u -i.o 




C 


>12 




NV 


3.6±1.4 


>2.3 


OI 


>6.8 




Sil 


<0.37 


<0.65 


Sill 


13±11 


>1.5 


Sill* 


1.7±0.7 


<0.71 


SilV 


>2.6 


>0.92 


A1II 


>0.16 
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<1000 cm -3 



we suggest that absorption feature A might be due to gas 
that has cooled and accumulated at the interface between a 
superwind bubble and the ambient IGM. 

In §3 we derived a lower limit to the mass of absorber 
A of > 10 s Mq. If we were to adopt an outflow velocity of 
~700 km s _1 (the typical blueshift of absorber A relative to 
the EELR: §3), the kinetic energy of absorber A would then 
be £ 5 x 10 56 erg. Could supernovae (SNe) inject sufficient 
energy into the ISM to power an outflow with this kinetic 
energy? Adopting a standard initial mass function, SNe yield 
~10 49 erg Mq 1 (B enson et al. 2003), meaning that a kinetic 
energy of ^ 5 x 10 56 erg would thus require a stellar mass of 
at least 5 x 1O 7 M0. This is compatible with the upper limit 
to the stellar mass of 4xl0 11 M Q for 2025-218 (Seymour et 
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al. 2007). Thus, if the absorber associated with 2025-218 is 
in outflow, it is energetically feasible that it is powered by 
SNe. 

The star formation rate would represent a further con- 
sistency check, because the rate at which SNe occur is also an 
important for powering an outflow. Unfortunately, the star 
formation rate in the host galaxy of MRC 2025-218 is un- 
known. Indeed, as yet there is no evidence for the presence 
of young stars therein (e.g. Cimatti et al. 1994). Sub-mm 
observations (e.g. Reuland et al. 2004) would be useful to 
estimate the rate of star formation in MRC 2025-218. 

Would the radio jets of MRC 2025-218 also be able to 
provide enough energy to power this outflow? Following Nes- 
vadba et al. (2006), the radio luminosity of MRC 2025-218 
(e.g. Carilli et al. 1997) implies a kinetic luminosity Lkin on 
the order of ~10 46 erg s _1 . The age of the radio source can 
be estimated by dividing its projected radius by its advance 
speed. For a projected radio source radius of 21 kpc (Car- 
illi et al. 1997), and assuming an advance speed of ~0.1c 
(e.g. Humphrey et al. 2007), we obtain an age of ~10 5 yr. 
Multiplying the kinetic energy of the radio source by its age 
then yields an energy of ~10 58 erg. Thus, it appears that 
the radio source is sufficiently energetic to have powered the 
gaseous outflow. However, we are not convinced that the 
anisotropic radio jets would be able to produce an outflow 
with the observed kinematic uniformity across such a large 
spatial extent. 

We note that in MRC 2025-218, we do not find evidence 
to suggest the extended ionized gas is in outflow. If absorber 
A is outflowing, it seems strange that the ionized gas at 
smaller radii does not show disturbed kinematics (Nesvadba 
et al. 2006) or morphological signatures suggestive of outflow 
(Zirm et al. 2005). 

4.2 Infall 

There has been much interest in recent years in the likely 
spectral signatures of the cosmological infall of gas onto 
galaxies at high-z (e.g. Haiman, Spaans & Quartaert 2000). 
Indeed, a number of recent investigations have revealed ev- 
idence for gas infall onto high-z active galaxies (Barkana 
& Loeb 2003; Weidinger, M0ller & Fynbo 2004; Humphrey 
et al. 2007). Different streams of infalling matter are ex- 
pected to collide at supersonic velocities, resulting in a quasi- 
spherical shock wave at or beyond the virial radius (e.g. 
Keshet et al. 2003). 

According to the model of Barkana & Loeb (2003), the 
accretion shock is expected to result in (i) an observable 
absorption feature in Lya, with a sharp edge marking the red 
side of this feature which corresponds to the infall velocity of 
Hydrogen crossing the shock, and (ii) a second Lya emission 
peak a few hundred km s _1 blueward of the absorption edge, 
due to the decreasing gas density towards larger radii and 
smaller infall velocities. Both of these features are present 
in the Lya profile of MRC 2025-218 (§3.1). We suggest that 
the radiation field produced by the post shock, cooling gas 
(e.g. DS96) might be responsible for ionizing the gas that 
has not yet crossed the shock (i.e., absorber A: see §3.5). 

Due to the relatively large expected radius of the accre- 
tion shock (~100 kpc: Barkana & Loeb 2003), we expect the 
absorption feature to have a large spatial extent and little 
variation in redshift across its full spatial extent. Absorber A 



associated with MRC 2025-218 is consistent with the infall 
model of Barkana & Loeb in this respect also. 

At the redshift of 2025-218, their model predicts a rate 
of gas accretion onto the host galaxy of ~200 Mq yr -1 . 
It is interesting to note that this accretion rate is in rea- 
sonable agreement with the order of magnitude infall rates 
that have been calculated for the EELR of (other) radio 
galaxies at z~2.5 (a few hundred M yr -1 : Humphrey et 
al. 2007). Analytical infall models imply a volume averaged 
Hydrogen density of pu ~20 times the cosmic mean den- 
sity, i.e., ~ 10~ 4 cm~ 3 , at the radius of the accretion shock 
(Barkana 2004). A volume filling factor of 10~ 7 would thus 
be required to reconcile the measured gas number density 
n e > 300 cm of absorber A with the predicted pu ■ 

Determining the metallicity of the absorbing gas is 
important for discriminating between infall and outflow. 
Naively, low metallicities would be expected if the absorbing 
gas is in infall (i.e., ^1 per cent solar: e.g. Shaye et al. 2003), 
while outflowing gas would be expected to have a metallic- 
ity of JjlO per cent solar. Unfortunately, in the case of MRC 
2025-218 our spectra have not allowed us to place any useful 
contraints on the metallicity of the absorbing gas. 

Using the analytical model of Barkana (2004) for the 
infall of gas onto a dark matter halo, we can estimate the 
column of H within a particular range of line of sight velocity 
(LOSV). The Lya absorber in 2025-218 is detected across a 
range in LOSV of ~270 km s" 1 (Fig 2: 4398.5 A to 4402.5 A), 
and thus we will consider only the H column from the LOSV 
of the accretion shock to a LOSV ^270 km s _1 blueward of 
this. We assume that the accretion shock radius is ~100 kpc 
and that the infalling gas has a velocity of 550 km s _1 when 
it reaches the shock (Barkana & Loeb 2003) . Based on these 
assumptions, we estimate that the infall rate of 200 Mq yr~ 
predicted from the model of Barkana & Loeb would imply 
a H column of roughly 3xl0 19 cm~ 2 within the relevant 
LOSV range. This estimate is likely to be a lower limit, 
because the post-shock cooling gas has not been considered. 
The H column density we have estimated for absorber A is 
not inconsistent with this value. 

Since the true systemic velocity is also unknown in other 
high-z radio galaxies, we suggest that, even when blueshifted 
relative to the line emission, the absorbers associated with 
some galaxies might in fact be in infall, rather than in out- 
flow as is usually supposed. Indeed, the HI columns of a 
number of the absorbers (e.g. Wilman et al. 2004) are in 
fair agreement with the HI column expected from the infall 
scenario of Barkana & Loeb (2003). By the same token, we 
suggest that the Lya absorbers sometimes associated with 
non-active galaxies at high-z (e.g. LAB-2: Wilman et al. 
2005) might also be related to the accretion of gas by the 
host galaxy. 



5 ABSORBER B 

The precise nature of absorber B is not clear. Its relatively 
large blushift with respect to the EELR emission of MRC 
025-218 allows us to reject infall. One possibility is that ab- 
sorption component B is in outflow from MRC 2025-218. 
A starburst driven super wind or a radio jet driven outlow 
would be energetically feasible (see §4.1). There is some ev- 
idence to suggest that MRC 2025-218 is located in galaxy 
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group (McCarthy, Persson & West 1992); we consider ab- 
sorption arising in an intervening starburst galaxy to be a 
possible alternative to outflow. The detection of NV in ab- 
sorption implies that young stars cannot be the only source 
of ionization for absorber B. Unfortunately, our measure- 
ments do not allow us to place any further constraints on 
the source of ionization. 



6 CONCLUDING REMARKS 

There is an extensive body of literature discussing absorp- 
tion lines associated with active galaxies (see e.g. Hamann & 
Ferland 1999; Crenshaw et al. 1999; and references therein). 
In essentially all of these investigations, the source against 
which the absorbing gas is seen is spatially unresolved, i.e., 
it has a size equivalent to that of the inner accretion disk. 
The overarching difficulty of these studies is to arrive at a 
firm conclusion about the nature of the absorber, i.e., it is 
not usually clear whether the absorber is a tiny plume of gas 
floating above the accretion disc, or a filament in the BLR 
located 1 pc from the accretion disc, or part of a galactic- 
scale super-wind, etc. Moreover, when discussing a spatially 
unresolved absorber, it is not immediately clear what such 
a small region of gas is representative of. 

By contrast, the study of spatially extended absorbers 
around galaxies is a relatively new approach to obtaining 
information about galaxies and their environments (e.g. van 
Ojik et al. 1997; Binette et al. 2006; Wilman et al. 2005; 
Wilman, Edge & Swinbank 2006). Giant nebulae of Lya 
emission (e.g. Villar-Martm et al. 2003) can provide a rel- 
atively bright and spatially extended background source by 
which intervening cold and warm gas clouds can be studied 
(van Ojik et al. 1997; Jarvis et al. 2003). If the absorber is 
found to be spatially extended, then its location in relation 
to the galaxy nucleus can be more readily ascertained. Un- 
fortunately, HzRG typically have relatively faint rest-frame 
UV continuum emission, and this normally precludes the 
detection of the multitude of other strong absorption lines 
that are expected to be present at wavelengths where line 
emission is very weak or absent. With only measurements of 
the Lya absorption line, and sometimes also CIV, it can be 
difficult to obtain definitive conclusions about the physical 
conditions and metallicity of the absorbing gas (see Binette 
et al. 2006 and references therein). We are of the opinion that 
galaxies that have both spatially extended Lya emission and 
very bright continuum emission hold much potential for fur- 
ther elucidating the true nature of associated absorbers. 

In this paper, we have presented a study of the ab- 
sorption lines in the rest-frame UV spectrum of one such 
galaxy: the z=2.63 radio galaxy MRC 2025-218. Thanks to 
the bright and spatially extended Lya emission, we have 
been able to trace Lya absorption across ~ 40 x 30 kpc 2 , 
where it shows remarkably little spatial variation in its prop- 
erties. This absorber is kinematically detached from the ex- 
tended Lya emitting gas, and has a covering factor of ~1. 
Based on these properties, we have concluded that the ab- 
sorber is outside of the EELR, i.e., at a distance of <;25 kpc 
from the galactic nucleus. 

In addition, MRC 2025-218 shows bright continuum 
emission from the active nucleus which has allowed us to 
measure metal absorption lines from a variety of atoms/ions 



and ionization levels. Although the continuum source, and 
hence the absorption lines, are not spatially resolved, we as- 
sociate them with the spatially extended Lya absorber since 
their redshifts and velocity widths are consistent, with the 
formal uncertainties. We find that the column density ratios 
are most naturally explained with photoionization by a hard 
continuum; an ionization parameter U~0. 0005-0. 005 is im- 
plied. We reject shocks and also photoionization by young 
stars, since neither is able to reproduce satisfactorily the 
measured column ratios. We have derived a lower limit of 
>10 cm -3 for the electron density of the absorber. The data 
do not allow useful constraints to be placed on the metallic- 
ity of the absorber. 

We have considered two possibilities for the nature of 
this absorber: the cosmological infall of gas, and an outflow 
driven by supernovae or the radio-jets. We find good agree- 
ment between the observed properties of the HI absorber 
and the properties of the HI absorption expected from the 
cosmological infall model of Barkana & Loeb (2003). It is 
also plausible that the absorber around MRC 2025-218 is in 
outflow. Discriminating between infall and outflow would re- 
quire knowledge of the systemic velocity of the host galaxy, 
and one way to achieve this could be to measure CI and CO 
lines using e.g. ALMA. A possible alternative could be to 
estimate the systemic velocity from the broad Ha emission, 
although a rather long exposure time would be required in 
order to reach a suitably level of high signal to noise (c.f. 
Humphrey et al. 2008). 
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